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^T) , Abstract 

(•~^ . We consider the group M of all polynomial v x v matrices U{uj) = ^I^qU^lo'', 

X"^ \ K = 0, 1, 2, ..., that satisfy equation U{uj)-D-U{uj)* = D with diagonal i^xv matrix D = 

^ I diag{— 1, 1, 1, ...1}. Here v > 2, U{uj)* = Yl'^=o ^*'^'' ^^^ U* is the Hermitean conjugate 

Q ■ of Ul. We show that the subgroup A^o of those U{uj) G Ai, that are normalized by 

^ , the condition C/(0) = /, is the free product of certain groups A4z where z is a z/-vector 

C^ . drawn from the set E = {z = columnjCi, ^2, •••G} I Ci = ^,z*Dz = 0}. Matrices in 

each Mz are explicitly and uniquely parametrized in the paper. Thus, every i/ x i/ 
matrix U = U{uj) G A^q can be represented in the form U = G^^ ■ Gz^ • ••• • Gz^ with 
v X v polynomial matrix multiples Gz, G -M-z,, z^ S H, for i = 1, ...ry, so that z^ ^ z^+i 
rS I for L = 1, ...7] — l and this representation is unique. The uniqueness includes the number 

C^ ' of multiples r] where 7/ = 0, 1, 2, ..., their particular sequence Gzi, Gz2, ■■■Gz^-, and the 

multiples themselves with their respective parametrizations; all these items can be 
defined in only one way once the U is given. 
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1. Let 

U{u;)=Uo + Uiu; + U2Uj^ + ... + U^uj'' (1) 

be a polynomial matrix (p.m.) of size z/i x 1/2- Here U^ for l = 0,1, ...k, are complex ui x 1/2 
matrices, cj is a scalar variable, and k = 0, 1, .... Under assumption that variable u takes on 
only real values, u* = u, we extend the Hermitian conjugatation operation Q to a p.m. U{uj) 
as follows: 

For i^ > 2, let D be the diagonal z/ x i/ matrix of the form D = diag{ — 1, 1, ...1}. In this paper 
we shall describe the group TW of all p.m.'s U^uj) of size u x u which satisfy the equation 

U{uj) ■ D ■ U{ujy = D, (2) 

i.e., the group of p.m.'s unitary in the indefinite metric of index 1. The description is effective 
in that every matrix in M is parametrized in a unique canonical form. The main results are 
stated in §§2, 3. This problem relates to the problem |jl]] of finding a factorization of a p.m. 

A{uj) = A{uj)* in the form 

A{uj) = U{uj) ■ C ■ U{u)* (3) 

with a p.m. U{u!) and a constant matrix C = C* , both of size v x u . Note that matrix C 
in (^ is not necessarily positive definite. As a special case, one might seek to factorize a 
constant matrix A{ijj) = C = C* , det C 7^ 0. In this case, without loss of generality, we can 
assume C to be a diagonal matrix 

C = diag{±l,±l,...,±l}. (4) 

The structure of the set of the solutions U{uj) of the equation (j^), where A{uj) = C , depends 
on how many +l's and — I's are on the diagonal of C. Indeed, if the elements of the diagonal 
of C are all of the same sign, then equation (|]), where A{uj) = C , coupled with the normal- 
izing condition t/(0) = I, has only the trivial solution U{uj) = /. The set of the solutions 
becomes non-trivial when C has both +l's and —I's on the diagonal. In this paper we 
consider the case of exactly one —1 with the rest of diagonal elements being +1. The same 
set of the solutions U will result if D is substituted with —D, a diagonal matrix that has 
exactly one +1 with the rest of diagonal elements being —1. These two cases exhaust the set 
of all indefinite matrices C of the form (^) in dimensions v <?>. Without loss of generality, 
we assume that the diagonal matrix C has the single —1 being at the top of the diagonal 

^For a constant vi x V2 matrix B, its Hermitean conjugate B* is defined as the V2 x vi transpose of the 
matrix whose elements are complex conjugates of the elements of B. If C/k 7^ in ([I]) then k is the degree 
of U{lo), denoted n = degC/(w). The following notations and conventions are assumed. Greek lower case 
denotes scalars, Latin lower case denotes column vectors of height v^ e.g., z — columnJCi, ^2, •■•Ci^}- Latin 
capitals denote matrices, i = V— T, / is the identity v x v matrix. The symbol 9 is read "such that" . The 
notation {Fi | F2} is used for the set with elements Fi satisfying the defining property F2. 



and the rest of the diagonal consists of +l's, i.e., C = D ioT the matrix D defined above. 
The task of finding a symmetric factorization (|^) recurs in many contexts, for example in 
the synthesis of linear optimal control in differential games . 

2. Introduce the sets: 
of vectors 

E = {z = column{Ci, ..., (^} | (^ = 1, z*Dz = 0} , 

A° = {d I d*z = d*Dz = 0}, for any z G H , 
of vector polynomials 

^z = {g = 9i^ + 92^"^ + ••• I 9i e A°}, for any z G S , 
and of scalar polynomials 

$ ^^ {0 = 0(^) I 0(0) = 0, <P{ujy = -0(u;)} . 
Then: 

1) the set of generatrices of the cone 

A = {x I x*Dx = 0} 
is parametrized in a one-to-one way by H ; 

2) A^ 7^ only if z/ > 2 and if ^f G A^ then g = Dg; 

3) an element G $ is of the form 0(c(j) = i >^ Pi.^\ where i = V— 1, the p^ are real 



l<t<K 



numbers, and k = 1, 2, .... 



For z ^S, (f) E ^, g E Az define the following u x u polynomial matrix 

G,(0, g) '=' D[z{<p - (1/2) ■ g*g)z* + {zg* - gz*)] + 1 . (5) 

The significance of the matrix in (^ is that, as will be shown, any matrix U E M. can 
be decomposed into p.m.'s of this form and a constant matrix. 

Define the sets of matrices: 

Mz = {G^(0, g)\(pe^,geA,} for any z G S , 



Mo = {f/M G M I f/(0) = /} , 

Af= {V eM\ degV^ = 0} , 

T''^{WeAr\W = diag{l, L}, Le{L\L*L = /(,_i),(,_i)}} . 

Clearly, every p.m. U{uj) G A4 may be represented uniquely in the form U{uj) = Uo{uj)-V 
where Uo{uj) G A^o, V eAf. 

The main result on the structure of Ai is the following 

Theorem 1. 1). For every 2; G S , the set Aiz is a group under matrix multiplication. 
The group TVIq is the free product ^ of groups M.z for all 2; G H. This means that any 
U = U{uj) G TWo can be decomposed in the form 

U = G,^{(l)i,gi)- ...■Gz^{(j)r,,grj) , (6) 

where Gz, G Aiz, for l = 1, ...rj. Moreover, if by aggregating the consecutive multiples in the 
sequence in (^ that belong to the same group we make sure that z,^ 7^ 2^+1, t = 1, ...rj — 1, 
then the obtained decomposition (|^) becomes unique for a given U . The uniqueness in- 
cludes the number of multiples rj, the particular sequence Zi ,Z2 , ...z^i in (^) with the sequence 
■Mzi , A^22 ; •••'^2,, of the corresponding groups, and the multiples Gz^ G Aiz, themselves, each 
of which is a unique p.m. for a given U{uj) G A^o- 

2). Multiplication in each Aiz, 2; G S, satisfies the condition 

Gz{^, g) ■ Gz{i^, h) = Gz{4>\^\ (1/2) • {h*g - g*h),g + h) , (7) 

where (f),ilj E ^, g,h E Az. 

3). Different groups M.z ore isomorphic; the connecting isomorphism is as follows 

W-Gz{<P,g)-W~^=GwMWg), (8) 

where ly G T, 2; G S, G $, (? G A^, and necessarily W~^ G T, Wz G S and Wg G Awz- 

4). The mapping {z,(f),g) \-^ Gz{(t),g) is one-to-one on the set 

{{z,(l),g)\ ZEE, G $, geAz, (p ^ 01 g ^ 0} 
IfU{uj) is decomposed in the form @ with z^ 7^ -Zt+i, '^ = 1, ■■■V — 1, then 
<legU{uj) = degG^,(0i,5fi) + ... + degG^^(0^,5f^J = deg(0i - glgi) + ... + deg{(l)r, - g*gr,) ■ 
Here z, G S, 0, G $, g, G A^^, i = 1, ...i]. 



5). The center ofAi^, z eE, is the set \l/j, = {Gz{4>,0) | G $} If u > 2, then the corn- 
mutant of M.^ coincides with the center of M.^- If^ = '^y then the group J^^ is commutative 
and coincides with \1/^ . 

3. Real cases. Let Ai' be the subgroup of Ai which consists of p.m.'s with real coef- 
ficients. Let S', $', A'^, z E E' he the real analogs of the sets S, $, A^ introduced in §2. 
Clearly, $' = {0}. As in §2 , let us define A^'^ and TWo- The structure of the group A^o is 
described in the following 

Theorem 2. If v = 2, then the group M.q is trivial. If u > 2, then the group A^o is 
not trivial and is the free product of groups A4'z for all z G S'. Groups Ai'z , ;z G S', are 
commutative. 

It is also interesting to describe the subgroup of p.m.'s in M. with real coefficients not 

with respect to variable u, but with respect to variable A = ia;. Note, that whereas uj* = u, 
we have A* = —A. Thus, a p.m. (|ID has real coefficients with respect to variable u if all U^ are 
real matrices; a p.m. (|I|) has real coefficients with respect to variable A if U^ are real for even t 
and if/i are real for odd l. As above, we introduce $" = {</) = 0(A) | 0(0) = 0, 0* = —</>}. The 
general form of such a is 0(A) = piA + p2A^ + psA^ + ..., where p^ are real numbers, and the 
sum is finite. In the u representation the general form of a G $" is i{piuj-\-p2Uj^+p3UJ^ + ...), 
with real p^. We define A^, ^A'., for z G S" and also A4"q. 

Theorem 3. 1). The group AI'q is a free product of groups Ai'l for all 2; G H" . 
2). If u = 2, then H" consists of two elements S" = {zi, Z2}, Z\ = column {1,1} and 
Z2 = column {1, —1}, so that the free product just mentioned contains two groups Ai], only, 
and each of these two groups is commutative; 

3). //z/ > 2, then each group Ai"^ is not commutative. The commutant of Ai'^ coincides with 
the center and is equal to {Gz{4>, 0) | G $"}. 

As an example of application of Theorem 3, let us present a general form of a 2 x 2 
p.m. U{X) G AIq of degree at most 2. We have: either U{X) = Gzj^{aX) ■ Gz2{PX) or 
f/(A) = Gz2{o:X) ■ Gzj^{PX). Here zi and Z2 are the two elements in E' when p = 2; aX and 
/5A are A-monomials in $" with real coefficients a and (3. Each monomial is of degree 1, if 
the coefficient is non-zero. Otherwise the monomial is zero. We can go further and represent 
these two types of U{X) in the per-component form. In the ffist case, we calculate that 



f/(A) = 
and in the second case, 
U{X)- 



-2a(3X'^-{a + (3)X + l 2apX'^ - {a - P)X 

2apX'^ + {a-p)X -2aPX'^ + {a + p)X + I 



-2af3X'^-{a + (3)X + l -2af3X^ + {a - f3)X 
-2apX^ -{a - I3)X -2apX^ + {a + I3)X + I 



Proofs. Let a\\h denote the existence of a linear dependence between vectors a and 6; we 
shall write a|||6 to mean that a\\b and, moreover, that the coefficients of the nontrivial linear 
combination of a and h can be chosen real. The proofs of Lemmas 1 to 2 below present no 
difficulties. 

Lemma 1. (a, 6 G A) =^ {{a* Dh = Q) ^^ {a\\h)) . 

Lemma 2. {ah* = ha*) <^=^ (a|||6) . 

Lemma 3. {XDX* = X*DX = 0) ^ (3y, zeE,a3 X = aDyx*) . 

The lemma is proved by applying Lemma 1 to the rows and, separately, columns of X. 

Lemma 4. Let y,z eE and a) XDz = 0; b) Wk e A° {XDk\\Dy). Then 3s G A°, r e 
A°, a3 X = D{rz* - ys* + ayz*) . 

Lemma 5. Let z^ ^S , l = 1, ...,1]. Then 

{DZiZl ■ DZ2Z* ■ ... ■ DZrjZ* = 0) ^^ (3io(l < io < ^) 3 z,Q = z,o+i) ■ 

The product of diadic matrices Dz^z*, l = 1,..."^, in the lemma is equal to {Dziz*)C,, 
where ^ = {zlDz2) ■ ... ■ (2;*_i-D-2^). But Dzi^;* ^ and the condition ^ = is equivalent to 
the existence of an 6, 1 < 6 < 77 — 1, such that zlDz^+i = 0. Since z^, 2:1+1 G A, it follows from 
Lemma 1 that 2J 2:1+1; hence, in view of the normalization, z^ = 2:^+1, and =^ implication 
is true. Obviously the <^^ implication is also true. 

Lemma 6. Let G $, g & A^, z Er. and suppose and g are not both zeros. 
Then 1) the leading coefficient of a p.m. Gz{4>,g) is proportional to Dzz*; 2) degGz{4>,g) = 
deg {(f) - g*g) > 0. 

The most laborous is the proof of 

Lemma 7. If degU{uj) > 0, then the degree of the p.m. U{uj) G M. may be decreased by 
a right or left multiplication by a p.m. of the form (^ (and real ifU{uj) is real). 

Proof of Lemma 7. Let U{uj) = Y1'^=o-^k.-l^^ , where X^,l = 0,1,...k, are complex 
(or real, or A-real, depending on the case) constant u x u matrices, and Xq 7^ 0, so that 
degU{uj) = K. It follows from (|^) that U{u)* ■ D ■ U{uj) = D. Extending definition of X^ for 
i > K to be null matrices, we obtain a family of equalities 

XoDx; + XrDx;_^ + ... + x^Dx; = {] {%) 

X* DX^ + XI D X^_i + ... + X; /} Xo = (10^) 

Here 7 = 0, 1, ..., 2/i; — 1. Having (Qq) and (IOq) we can apply Lemma 3 to X = Xq. For the 
leading coefficient Xq, this yields its diadic representation Xq = aQDyz* for some vectors 



y,z & E , and a number ao ■ Note that ao 7^ because deg U{uj) = k . The y and z will 
be our candidates for the subscript of the p.m. of the form (H) which should decrease the 
degree of Uiu) after multiplying U{uj) on the left or on the right, respectively. Specifically, 
we will prove the lemma if we show that at least one of the two possibilities holds: 

I) 30 e $, 3(7 G A, 9 deg [U{uj) ■ G^icj), g)] < k, 

II) 3^^ e<^,3heAy3 deg [Gyiil^, h) ■ U{u)] < k . 

(In these statements, sets S, $, A^;, A^, should be appropriately modified in the cases of reals.) 

For integer positive numbers r, /i, and ^, such that t < fi and t < C,, consider the follow- 
ing conditions: 

A) 3a(, ,0<6<r — l,3Xt = a^Dyz* , i.e., the first r leading coefficients of p.m. Uiu) are 
diadic matrices proportional to Xq\ 

B) Xi^Dz = for all i, r < 6 < /i — 1, and X^Dz 7^ , i.e., the fi — t coefficients of p.m. 
U{uj), that follow the last coefficient Xr-i mentioned in A), turn into 0, when multiplied on 
the right by the vector-column Dz, but this does not haappen for the /i — r + 1st coefficient; 

C) y*Xi_ = 0, for all t, r < t < ^ — 1, and y*X^ 7^ , i.e., the ^ — r coefficients of p.m. U{uj), 
that follow the last coefficient Xr-i mentioned in A), turn into 0, when multiplied on the 
left by the vector-row y*, but this does not happen for the ^ — r + 1st coefficient; 

D) 3 s e A° , r G A° , a^ 9 X^ = D{rz* — ys* + a^-yz*) , i.e., the coefficient X^- is the sum of 
three diadic matrices as stated. 

Consider the largest r 3 A). Since Xq = aoDyz*, the r is at least 1. It can not be larger 
than K, though. Hence, 1 < r < k. Let yU 9 B), ,^ 9 C). Clearly, t < fi,C, < n- The way the 
proof proceeds further depends of whether the two inequalities 

/i>2rand^>2r (11) 

both hold or not. First we consider the easier 



Case 1: at least one inequality in (|TI]) fails. Suppose, for example, that /i < 2r. We will 
then verify I). Equality (9^) and Lemma 2 imply that there exists a real number p such that 

aoDy + ipXf,Dz = 0. (12) 

Using (|12D it is easy to verify that deg [U{uj) ■ {Dz{ipuj^)z* + 1)] < k, i.e., as stipulated in I), 
the degree decreases when U{uj) is multiplied on the right by Gz{(p,g) with (p{uj) and g^uj) 
taken here as (p{uj) = ipuj^ , g{u) = 0. (When z, y and X^ are real, since p is real, (0) implies 
that iao is real and hence all elements of iXo are real. It follows, that if U{uj) has only real 
coefficients, inequality p < 2t can not occur. In case of reals with variable A = iu, inequality 
p < 2t can only occur for odd p.) Analogously we verify II) if the second inequality in ([Til) 
fails. 



Case 2: both inequalities in ([TT|) hold. This case will be more laborous. First, we will 
verify D) by using Lemma 4 where X = Xr . The conditions a) and b) in Lemma 4 are obvi- 
ously satisfied. Let us verify the condition c). Pick a vector A; G A° and denote v = XrDk. 



Condition c) is then rewritten as v\\Dy. To verify the latter, we are going to use Lemma 1 
with a = Dy, b = v. Obviously a G A, so we have to check that also I') 6 G A, i.e., v*Dv = 0, 
and that IF) a*Db = 0, i.e., y*v = 0. 

Multiply the equality (92r) by Dk on the right and by k*D on the left. It follows from 
(|TT|) that on the left-hand side all the summands but one will turn to zero. In other words, 
r) holds. We also have y*v = {y*Xr)Dk = 0, that is, 11') holds. 

Now, as the pre-conditions of Lemma 4 are satisfied, we obtain representation 

Xr = D{rz* - ys* + ayz*) . (13) 

Observe, that vectors s and r can not be both zero in (|1^), because it would have contra- 
dicted to the definition of r which could have been possible to increase in such a case. Let 
us assume that s 7^ in (|13[) and show I). If r 7^ in (|1^) we can similarly show II). 

Define w = X2rDz,p = a^Dy. Using Lemma 1, we will show that w\\p. Obviously 
p*Dp = and we will also show that I") w*Dw = II") p*Dw = which are the precondi- 
tions in Lemma 1 applied to a = w, b = p. 

Using (lO^r) immediately yields I"). To obtain II"), we begin with (102r) which we 
multiply on the right at Dz. This yields aQzy*w = which implies 11"). Now we can apply 
Lemma 1 and obtain 

w = (ao + ipo) -p, (14) 

where o"o and po some real numbers. (In the case of reals with respect to variable u the 
resulting po will be zero, in the case of reals with respect to variable A the resulting po will 
be zero for ?????). 

We can now use the obtained expressions for substituting in (92r) 
Xo = pz*, X2r = (o-Q + ipo) -p X^ = -pso + nz* 

where we denoted n = D{r + a^y) and Sq = s/a^. We then obtain an equation pp*{2ao + 
SqDsq) = from which, taking into account equality SqDsq = |soP, we obtain o"o = 

-(l/2)kop. 

We are now ready to verify I). Set (piuj) = ipu'^^ and g{uj) = du^, where rf G A° and 
real p are constants to be defined. The definition should satisfy condition 

Xo + X2rDzz*{ip - (l/2)|rf|2) + XrD{zd* - dz*) = pz*e = 0. (15) 

Parameter e which is defined in (|TB|) is equal to 

e = l + sld+ (ipo - (l/2)|son ■ (ip - (l/2)Mn. 
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Our goal is to set parameters p and d so that e will become zero. We are looking for d of the 
form d = 9so, where Oq is a constant to be defined. It follows, that e = ei — i ■ (l/2)|soP ■ 62. 
where ei = |1 + (1/2)6' ■ |soPP — P • Po and €2 = po ■ |6'p — 2 ■ Im6' + p. If po = 0, we set p = 
and 9 = — 2|so|~^ to obtain e = 0. If Po 7^ 0, then by setting p = |1 + {1/2)6 ■ |soPP/po we 
obtain ei = 0. Such setting of p results in 

62 = Po'{l + \0\'- [pI + ((l/2)|^or)1 - 2[po ■ Im^ - iil/2)\s\') ■ ReO]} (16) 

Now setting 

e=(ipo-(l/2)|soP)-(p^+((l/2)|aoP)^)-^ 

results in €2 = as required. Lemma 5 is proved. 

Proof of Theorem 1. Let us first establish (|^). We have 

G,(0, g) ■ G,{^, h) = {DA + I){DB + I) = DADB + D{A + B) + I, 

where A = z{(f)—{l/2)-g*g)z* + {zg*—gz*) and B = z{ii — {1/2) ■h*h)z* + {zh* — hz*) according 
to (^). Using identities z*Dz = 0, z*Dh = 0, and g*Dz = 0, that follow from the definitions, 
and multiplying each summand in A by each summand in B on the right, we establish that 
the only non-zero term in the resulting sum for DADB is —Dzg*Dhz* = —Dz{g*h)z* where 
Dh is substituted with h since h E A^. Using the identity, 

{g + h)*{g + h)- {h*g - g*h) = g*g + h*h + 2g*h , 

the sought product then becomes 

G,(0, g) ■ G,{iIj, h) = -Dzg*hz* + D{A + B) + I = 

= D[z{<P + ^ - (1/2) ■ {g*g + h*h - 2g*h))z* + {z{g + h)* - {g + h)z*) 

= 6-^(0 + ^ + (1/2) • {h*g - g*h),g + h) 

as required in (0). Similar and even simpler arguments show that for U = Gz{(f), g) @ holds. 
Hence sets A^^ for 2; G S are indeed groups under the matrix multiplication and TVI^ C TVIq 
for any 2; G S. 

It is easy to see that deg Gz{(t>, g) = deg(0 — g*g) unless both (j) and g are identical zeros 
in which case Gz{4>,g) is identical to the matrix /. Hence if Gz{4>,g) is not identical to 
/, then the leading coefficient of Gz{4>,g) is proportional to the diadic matrix Dzz*. If a 
product of matrices Gz{4>,g) is formed, like that in the right-hand side of (^ with adjacent 
multiples belonging to different groups Aiz, then the leading term of this product will be 
proportional to the product of the diadic Dzz* as in the left-hand side of the statement in 
Lemma 5. According to that lemma, the leading term can not degenerate to zero. Because 



if it did, then a pair of adjacent multiples would have belonged to the same group Aiz- 
Therefore, the products like those in the right-hand side of (P) can never degenerate to the 
identity matrix /, given that the number of multiples is not zero and that adjacent multiples 
belong to different groups Ai^- Hence the product of group A^^ where z runs over set S, 
we temporarily denote this product as A^, is free (i.e., has no relations, see, e.g., [^ for the 
definition). Lemma 7 tells that A^o C Al and obviously we have Ai C Alo, thus Ai = Alo 
and statement 1) of theorem 1 is proved. 
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